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Physical transformations between quantum states 

Zejun Huang, Chi-Kwong Li, Edward Poon, Nung-Sing Sze 



^ ! Abstract 



Given two sets of quantum states {Ai, . . . , Af^} and {Bi, . . . , -Ba,}, represented as 
sets as density matrices, necessary and sufficient conditions are obtained for the ex- 
istence of a physical transformation T, represented as a trace-preserving completely 
positive map, such that T{Ai) = Bi ior i = 1, . . . ,k. General completely positive 
maps without the trace-preserving requirement, and unital completely positive maps 
' (— I 'i transforming the states are also considered. 
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1 Introduction and Notation 

(N. 

> ■ 1.1 Introduction 
l> 

ly-v ■ In quantum information science, quantum states with n physically measurable states are 

^i ' represented hy n x n density matrices, i.e., positive semidefinite matrices with trace one. In 

Q I particular, pure states are rank one density matrices, while mixed states have rank greater 

^ ! than one. We are interested in studying the conditions on two sets of quantum states 

{Ai, . . . , y4fc} and {-Bi, . . . , -B^} so that there is a physical transformation (a.k.a. quantum 

k> ■ operation or quantum channel) T such that T sends A^ to Bi for i = 1, . . . ,k. 

H ■ To set up the mathematical framework, let M^ „ be the set of m x n complex matrices, 

and use the abbreviation M„ for M„,^. Denote by x* and A* the conjugate transpose of 

vectors x and matrices A. Physical transformations sending quantum states (represented as 

density matrices) in M„ to quantum states in M^ are trace-preserving completely positive 

(TPCP) maps T : M„ — )■ M^ with an operator sum representation 

r 

t{x) = J2FjXf;, (1) 

where Fi, . . . ,Fr are m x n matrices satisfying Y7j=i ^j^j ~ ^n] see [3], [5], [TJ §8.2.3]. So, 
we are interested in studying the conditions for the existence of a TPCP map T of the form 
([ID with ^^^^^ F*Fj = In such that T{Ai) = Bi ior i = 1, . . . ,k. 



We also consider more general types of physical transformations (completely positive (CP) 
linear maps) without the trace-preserving assumption, i.e., not requiring Yl^j=i-^j^j — ^n- 
Such operations are also considered in the study of quantum information science; see [TJ 
§8.2.4]. Furthermore, in Section 4 we consider unital completely positive maps which are of 
interest in the theory of C*-algebras. Such CP maps are dual to the trace-preserving ones 
and send the identity matrix to the identity matrix, i.e., they satisfy Y7j=i ^j^j ~ ^m- 

In Section 2, we study physical transformations on qubit states, i.e., quantum states on 
M2. Section 3 concerns physical transformations sending general states to general states, 
and Section 4 concerns more general transformations acting on pure states. 

1.2 Notation 

We conclude this section by defining additional notation and recalling some terminology that 
will be used later. Given a matrix M (which we may alternatively denote as (Mij), to focus 
on its entries), we write M* for the transpose of M, and M for the matrix whose (i, j)-entry 
is the complex conjugate of M^j. The Hadamard product (or Schur product) of two m x n 
matrices A and B is the m x n matrix Ao B whose (i, j)-entry is given by AijBij. (So, the 
o symbol denotes entry-wise multiplication.) A correlation matrix is a positive semidefinite 
matrix with all diagonal entries equal to 1. 

Suppose a matrix A has the spectral decomposition A = XlfeLi ^k'^k^l for some orthonor- 
mal eigenvectors v^- One possible purification for A is the vector ^11^=1 \f^'^k®Vk'i the most 
general form for purifications of A are vectors of the form = 'Yll^=i^f^'^k®Wvk G C^^C**, 
where TV is a partial isometry from C"^ to C^. Note that, for any purification of A, the 
partial trace of 00* over the second system is precisely A, and one can actually take a more 
abstract point of view and define purifications to be those vectors possessing this property. 
(Recall that the partial trace of S ® C e Mm ® M^ over the second system is just -B(trC), 
and one extends linearly to define the partial trace on all of Mm ® Mr.) 

2 Qubit states 

In this section we focus solely on qubit states (2x2 density matrices). Recall that the trace 
norm || ■ ||i of a matrix X is the sum of its singular values. The following interesting result 
was proved in pLj; see also 0. 

Theorem 2.1. Let Ai,A2, -81,-62 G M2 be density matrices. There is a TPCP map sending 
Ai to Bi for i = 1,2 if and only if \\Ai — tA2\\i > \\Bi — ti?2||i for all t > 0. 

The proof in pj is quite long. In the following we give a short proof of the result, and give 
another condition that is much easier to check (condition (c) in Theorem 12. 2 p by making the 



following reduction: if rank Ai = 2, then we can find c > so that Ai = A\—cA2 is a positive 
semidefinite matrix of rank one. Then we simply replace Ai, Bx by Ai, Bi = Bi — CB2, since 
a TPCP map sending Ai to -B, exists if and only if there is a TPCP map sending Ai to Bi 
and A2 to i?2- We may then repeat the process by considering A2 = A2 — cAi. 

So, by taking linear combinations of Ai, A2 (and the corresponding combinations of Bi, 
B2), we may assume that Ai = Xix* and A2 = X2X2. We have the following. 

Theorem 2.2. Let Ai = xix*, A2 = X2X2, Bi, B2 G M2 be density matrices. The following 
conditions are equivalent. 

(a) There is a TPCP map sending Ai to Bi for i = 1,2. 

(b) ^il + ty-4:t\xlx2\'^ = \\Ai - Malli > \\Bi - t^alli for all t > 0. 

(c) \xlx2\ = llv^v^lli < llVB^v^lli. 

Note that condition (c) is of independent interest, for it relates the fidelity between the 
initial states with the fidelity ||-\/&i"-\/52||i between the final states Bi, B2, and can be 
generalized to give a necessary (but not sufficient) condition for the existence of a TPCP 
map sending k initial states to k final states (see equation IQ later, also [2]). 

Proof Note that for X e M2, ||X||f = tr (XX*) + 2| det(X)|. One can readily verify the 
first equality in (b) and the first equality in (c). 



(a) => (b). Suppose T is TPCP. If A = A^ 
inite, then 



A_ where A^ and A_ are positive semidef- 



lin^)lli < \\nA^)\\i + r(A_)|K = trT(A+) + trT(A 



tiA, +try4_ 



l^lli- 



Thus \\Bi -tB2\\i = \\T{Ai - tA2)\\i < \\Ai - tA2\\i for alH > 0. 

(b) =^ (c). Suppose one of the matrices Bi and B2 has rank 1. Without loss of generality, 
we may assume that B2 = 1/21/2 ■ By condition (b), for t > |/|i?i|/2, we have 

{1 + tf - U\xlx2\^ > \\Bi - ty2y*2\\l 

= tr ((fii - ti/21/2)') + 2| det(Si - ty2y*2)\ 
= t^ + 2t-4t{y*2Bm)+l 



for a constant 7 G M. Thus, |x];x2p < 1/2-^11/2 



^^v^" 



'2||l- 



Suppose both Bi and B2 are invertible. Choose t so that det(i?i) = det(t-B2)- Applying 
a suitable unitary similarity transform, we may assume that Bi — ti?2 is in diagonal form so 
that 



Bi 



hi c 
c I -hi 



B, 



62 c/t ' 
c/t 1-62 



Then 



det(5i + tB2) - I det(fii - t^a)] 

[(61 + tb2){l + t-h- th) - 4|c|2] - (61 - t62)((6i - th^) - (1 - t)) 



2{6i(l-6i) 



+ t\h{l-h)-\c\'/e)} 



2(det(5i) + det(t52)) 

4det(/B;v^) 
4tdet(v/5^\/5^)- 



Hence, 



because t satisfies det(i?i) = det(ti?2) 



4tdet(Vfi^V^) 



det(Si + tB2) - I det(5i - 182) 
By condition (b), we have 

(l+t)2-4t|x*X2p 
> tr((5i-tS2)2) + 2|det(5i-t52)| 
= tr ((5i + t^s)') - 2ttr (^1^2 + ^s^i) 

+2 det(5i + t52) - 2 det(5i + t^s) + 2| det(5i - ^^2) | 
= (tr (5i + t52))2 - 4ttr {B1B2) - 2det(Si + 182) + 2| det(5i - tB2)\ 

= (l+t)2-4t[tr(fii52) + 2det(/B^/B^)] by® 

= {l+tf-At\\^/B[^2\\l 

Thus, ||a/B7v^5^||i > |x^X2|, and condition (c) holds. 

(c) ^ (a). Note that ||X||i = max{|trXiy| : W is unitary}, so there exists a unitary 
V G M2 such that \ti \fWi\fB2V\ > |a;^a;2|. If we write \fBi = [yi\y2] and \/B2V = [zi\z2], 



and set y 



yi 

y2 



G C and z 



G C^, then this inequality implies that \y*z\ > |x^a;2|. 



Set (5 = 1 if y*z = 0; otherwise let 6 = {xlx2) / {y* z) . Then the 8 x 2 matrices 



X 



Xi X2 





and Y 



y ^-2 



satisfy X*X = Y*Y (note that yiy^ + y2yl = Bi and ziz^ + ^2-^2 = -^2, so taking the trace of 
these equations shows that y and z are unit vectors), so there exists a unitary U such that 
UX = Y. Regard the first two rows of U* as [F^F2F^F^]. Then the map 

X ^ FiXF* + ■■■ + FiXF* 

is the desired TPCP map. D 

Remark. Consider the problem of the existence of a TPCP map T such that T{Ai) = Bi 
for i = 1, . . . ,k, for given density matrices Ai, . . . ,Ak,Bi, . . . ,Bk G M2. Evidently, we can 



focus on the case when {Ai, . . . , Ak} is a hnearly independent set. li k = 1, then the map 
defined by T{X) = (trX)Si is a TPCP map satisfying the desired condition. Theorems 
12.11 and 12.21 provide conditions for the existence of the desired TPCP map when k = 2. 
If fc = 4, then {Ai, . . . , A4} is a basis for M2. There is a unique hnear map T satisfying 
T{Ai) = Bi ioT i = 1, . . . ,4. It is then easy to determine whether T is TPCP by considering 
its action on the standard basis {En, E12, E21, E22} for M2. One simply checks whether 
tr (Ell) = tr {E22) = 1, tr (-^12) = tr (-E21) = 0, and whether the Choi matrix 

T{En) T{E,2) 
T{E2i) T{E22) 

is positive semidefinite; see [3]. The remaining case is when fc = 3. Again, we can re- 
place Ai,A2, A3 by suitable linear combinations (and apply the same linear combinations to 
Bi,B2,B3 accordingly) and assume that Ai = XiX* for i = 1,2,3. We have the following 
result. 

Theorem 2.3. Suppose Ai = XiX*,Bi G M2 are density matrices for i = 1,2,3 such that 
Ai,A2,As are linearly independent. Let X3 = «ie**^Xi +026**^3^2 with ai,a2 > 0, ti,t2 ^ 
[0, 27i), and 



B. 



1 



-{B3-alBi-alB2). 



2aia2 

Then there is a TPCP map sending XiX* to Bi for i = 1,2,3 if and only if there exists 
C e M2 such that 



ti (CC*) = 1 + 1 det{C)\' < 2, tr ^JB2C^B^ = e'^'^-''>xlx2, and B^ = Re^/B2C^/Bi. 

(3) 

Proof. First, consider the forward implication. Note that T is a TPCP map sending XiX* to 
Bi for i = 1,2 if and only if \x\x2\ < \\^/B^^/B^\\l. li we write T(X) = J2'j=iFjXF* and 
FjXi = yij, note that Yi = [yn . . . yir] must equal y/B^W* for some isometry Wi G M^m- 
Writing Re A= {A + A*)/2, we have 

r r 

T{x3xl) = ^FjX^xlF* = ^ {a\yijy{j + aly2jy*2j + 2Reaia2e'^^^-^^^y2jylj) 

= alBi + alB2 + 2aia2Re e'^*^-*'^Y2Y* 
= aJBi + ^2^2 + 2aia2Re>/B^C >/bI, 

where C is a contraction and tr y/B^C\^B^ = e^^^'^~^^^x\x2- Note that C is a contraction 
if and only if the largest eigenvalue of CC* is bounded by 1, which is equivalent to the 
inequalities: 

tr {CC*) < 1 + det{CC*) = 1 + | det(C) ^ < 2. 



Suppose the first inequality is a strict inequality. Consider the subspace 

S = {X eM2 :Rev/5^Xv^ = 0.tr(v/5^Xv/fi^) = 0} C M2. 

Then we may replace C hj C + X with X E S so that ||C + X|| = 1, and the new solution 
C will satisfy the equality tr (CC*) = 1 + det(CC*). 

Conversely, suppose there exists C satisfying condition (3). Write y/Bi = [^11^12], 
\fB^C = [2/21 1/22], and \fB^^{I — CC*) = [2/232/24]- Then the inner product of the two 



unit vectors e**^xi and e**^X2 equals that of the unit vectors 
is a unitary U G Mg such that 

U 



2/11 




2/21 


2/12 




and 


2/22 

2/23 


_0_ 




.2/24. 



. Thus, there 









2/11 


2/21" 


Oe 


e**2x2 
Oe . 


= 


2/12 

O2 

.O2 


2/22 

2/23 
2/24. 



Let the first two rows of U* be [F* F; F* F*]. Then the map T{X) = J2.=i FjXF* satisfies 



T{A, 



Bj for i 



1,2,3. 



D 



Note that condition (|3]) can be verified with standard software. In fact, if we treat C as 
an unknown matrix with 4 complex variables (that is, 8 real variables), then the last two 
equations translate to 5 independent real linear equations. By elementary linear algebra, the 
solution has the form C = Cq + XiCi + X2C2 + x^C^ for 4 complex matrices Cq, Ci, 6*2, C*3 
in M2, and 3 real variables Xi^X2-,x^. Then we can substitute this expression into the 
first equation to see whether the first nonlinear equation (of degree two) is solvable. In 
fact, we can formulate the first equation as an inequality: tr(CC*) < 1 + |det(C)p < 2. 
Then standard computer optimization packages can decide whether there exist real numbers 
Xi,X2,X:i satisfying the inequalities. 



3 General states to general states 
3.1 Moving beyond qubits 



A natural question is whether or not Theorem 12. 21 can be generalized to non-qubit states, i.e. 
states on M„ where n > 2. The equivalence of (a) and (b) in Theorem 12.21 does not hold for 
density matrices with dimension greater than two (a counter-example may be found in [4]). 
On the other hand, it is known (see [21 Lemma 1], for example) that the equivalence of (a) 
and (c) holds for density matrices of any dimension — provided the initial states Ai, A2 are 



6 



pure, i.e. have rank one. (See the example below.) This illustrates two points. First, results 
for states of arbitrary dimension appear to be more readily attainable when the inputs are 
restricted to be pure. Second, this shows why the situation is easier for qubit states: for 
qubits, one can always perform the reduction described before Theorem 12.21 to reduce to the 
case where the input states are pure, whereas this cannot be done in general for non-qubit 
states. 

Example. Note that ||v^]47v^^||i < ||-\/57-\/fi2||i does not imply ||Ai— tA2||i > ||-Bi — ti?i||i 



for all t > if Ai and A2 are not of rank one. For example, let Ai 
diag (1/3, 2/3) and 



diag(4/5,l/5),A2 



Bi 



1/4 ^3/4 
73/4 3/4 



Bo 



Then 



while 



%VA2l 



0.8815 < 0.9659 



1/2 1/2 
1/2 1/2 



'BiJB^Wi 



\Ai - 5A2II1 = 4 < 4.1641 = \\Bi - 53; 



2 1- 



So, what more can be said if we impose the additional restriction that the initial states are 
pure? Well, if we also assume that the final states are pure, we have the following interesting 
result from [2| Theorem 7]. 

Theorem 3.1. Let Xi G C" and r/i G C'" be unit vectors for i = 1, . . . ,k. Let X = [xi\ . . .\xk\ 
and Y = [yi\ . . .\yk]. Then there exists a TPCP map T such that T{xiX*) = yiH*, i = 1, . . . ,k 
if and only if X*X = M o Y*Y for some correlation matrix M G M^. 

Note this gives a computationally efficient condition to check if the matrix Y*Y has 
no zero entries. We will use this result as a model to generalize in the rest of the paper, 
considering the most general situation ffist in the next subsection (where we obtain a result 
which allows us to derive the above theorem as a special case), and then, in the subsequent 
subsection, we consider keeping pure input states, but relax the condition that the final 
states be pure. The final section examines how this theorem changes when the maps are not 
necessarily trace-preserving. 



3.2 Mixed states to mixed states 

In this subsection we consider the difficult problem of characterizing TPCP maps sending 
k initial states to k final states (not necessarily of the same dimension), starting with the 
general case, and then considering special cases that are more tractable. The following 
theorem is rather technical, but it does provide a useful framework for the most general 



situation, and can be readily applied to quickly derive existing results under more specialized 
circumstances. The multiple equivalent conditions reflect various approaches and serve as a 
segue between different viewpoints and lines of attack on a problem. Note that we ignore 
zero eigenvalues when using the spectral decomposition in the theorem's statement so as to 
eliminate redundancies, thus preventing matrices from becoming artificially large. 



Theorem 3.2. Suppose Ai, . . . , A^ G M„ and Bi, . . . ,Bk G Mm are density matrices. Using 
the spectral decomposition, for each i = 1, . . . ,k, we may write Ai = XiDfX* and Bi = 
YiDfY* , where Xi,Yi are partial isometrics, and Di G Mr^, Di G Mg^ are diagonal matrices 
whose diagonal entries are given by the square roots of the positive eigenvalues of Ai, Bi 
respectively. The following conditions are equivalent. 

(a) There is a TPCP map T : M„ — )■ Mm such that T{Ai) = Bi for i = 1, . . . ,k. 



(b) For each i 



k and j 



, Ti, there are Si x s matrices Vij such that 






V--V* = I 



and the {p,q) entry of the r^ x Vj matrix {DiX*XjDj) equals tr {V*pD*Y*YjDjVjq) . 

,1' 



(c) There are vectors x,- 



Xli 






G (C")'' and vectors yji 



y 



j« 



y 



j«. 



G {C"")' fori = l,...,k 



and j = 1, . . . ,r such that Ai = X]j=i ^ji^*ji! Bi = X]j=i St=i y%i.y]d* '^^^ there is a unitary 
U G Mms satisfying 



U 



Xii ■ ■ ■ Xri Xi2 ■ ■ ■ Xrk 

Urns— n ' ' ' ^ms—n ^ms—n ' ' ' ^ms—i 



y\i 



yrk 



yii ' ' ' yrk 



Proof, (a) =^ (b). Let Cj denote the vector with 1 in the ith position and in the other 
positions. Note that AA* < BB* in the Loewner order (that is, BB* — AA* is positive 
semidefinite) if and only ii A = BC for some contraction C . As in equation ([T]), we may use 
the operator sum representation for a TPCP map to write T{Ai) = Yl^i=i Pi^iPi fo^ some 
m X n matrices Fi. Thus 



(r,A)(K,A)* = 5, = T{A,) = y^FiX,DJX*F: > {FiX,D,e,){FiX,D,e,y 



1=1 



YiDic\j for some vectors c'^ G 



for any i,j,l, whence FiXiDiCj 

Since T{Ai) = Bi it follows that J^^i VijWj = h^ 



Let Vii 



\c^-\ 



■ ■ \C1 



8 



The trace-preserving condition X]f=i ^i^i — ^n implies that there is a unitary matrix 
U G Mms whose first n columns are given by [Fj* . . . F*~\ . The rest of (b) follows by 
noting that the inner product of any two columns of the ms x (ri + • ■ ■ + r^) matrix 



X 



XiDi ... XkDk 

^m.^i — T). • • • ^m.et — T) 



must equal the inner product of the corresponding two columns of UX. 

(b) =^ (c). Let r = maxjrj. Set Xji = XiDiCj if j < r^ and Xji = otherwise. Let 
y'-j- = YiDiVijCi for / = 1, . . . , s if j < r,, and set ylj = otherwise. Then the summations 
to Ai and Bi are clearly satisfied. Finally, the last condition of (b) implies that the inner 
product of Xpi and Xgj equals the inner product of ypi and ygj, and this entails the existence 
of a unitary U satisfying the final condition of (c). 

(c) =^ (a). Let [-Fj* . . . F*~^ be the first n rows of U*, and define T by T{X) = 
^"j^^FjXF*. The resuh follows. D 

The conditions (b) and (c) are not easy to check. It would be interesting to find more 
explicit and computationally efficient conditions. Nonetheless, one can use the above theorem 
to deduce Corollary 10 in [2J for TPCP maps from general states to pure states. 

Corollary 3.3. Suppose Ai, . . . ,Ak G M„ and Bi = yy^, . . . ,Bk = yuvl ^ ^m are density 
matrices. For each i = l,...,k, write Ai = XiDfX* such that Di G Mr- are diagonal 
matrices with positive diagonal entries. The following conditions are equivalent. 

(a) There is a TPCP map T : Mn — )■ M^ such that T{Ai) = Bi for i = 1, . . . ,k. 

(b) For each i = 1, . . . ,k and j = 1, . . . , Tj, there are vectors Vij such that ^ 
and the {p,q) entry of the r^ x Vj matrix {DiX*XjDj) equals v*pVjqy*yj. 

(c) For each i = 1, . . . ,k and j = 1, . . . , Tj, there are vectors Vij such that Y17=i 
and a unitary U satisfying 



j=i^.j^^J 



VijVij 



U 



XiDi ■ ■ ■ XkDk 
... 



[vn^yi ■■■ vir^®yi ■■■ Vki®yk ■■■ Vkr^® Vk] ■ 



When all Ai and Bi are of rank one, the above result reduces to the following result. 

Corollary 3.4. Suppose Xi, . . . , x^ G C" and yi, . . . ,yk G C'" are unit vectors. The following 
conditions are equivalent. 

(a) There is a TPCP map T : M„ — )■ Mm such that T{xiX*) = yiy* for i = 1, . . . ,k. 

(b) There exist Dp = diag (f ip, . . . , Vkp) for p = 1, . . . ,s satisfying 

s s 

j=i p=i 



(c) There is a correlation matrix C G Mk such that 

{x*Xj) =Co{y*yj). 

(d) There are unit vectors vi, . . . ,Vk G C'^ and a unitary U G M^s such that 



U 



Xi 



Xk 



Urns— n ' ' ' ^ms—n 



[f 1 ® yi ■ ■ ■ VkiS) Vk] ■ 



Note the equivalence of conditions (a) and (c) above is just Theorem 13. II 

3.3 Pure states to mixed states 

Next, we turn to TPCP maps sending pure states to possibly mixed states, and give a number 
of necessary and sufficient conditions for their existence. This problem was also considered in 
[2] using the concept of multi-probabilistic transformations. We instead rely on purifications 
of mixed states, with the aim of generalizing Theorem 13.11 

Theorem 3.5. Suppose xi, . . . ,Xk G C" are unit vectors and Bi, . . . ,Bk G M^ are density 
matrices. Then there is a TPCP map T such that T{xiX*) = Bi for i = 1, . . . ,k if and 
only if there exist purifications yi of Bi such that X*X = Y*Y , where X = [xi| . . . \xk] and 
Y = [yi\... 



Proof Suppose there is a TPCP map T such that T{xiX*) = Bi. Write T{A) = ^^^^^ FjAF*. 
Since T is trace-preserving, [F* F2 ... F*~^ has orthonormal rows and can be extended 
to a unitary matrix U* G Mmr- Define yi = U{xi © Omr-n)- Write 



yii 

Uri 



G (C 



X 



Xi 



Xk 



(4) 



Then FjXi = yji, and Bi = T{xiX*) = J^^j^iVjiyji, so yi is a purification of Bi. Moreover 
Y*Y = {UXyUX = X*X = X*X as desired. 

Conversely, suppose we have purifications yi oi Bi, written as in (jl]) with Bi = Yl^j=i VjiVji- 
If y*y = X*X = X*X, then, since Y and X have the same dimensions, there exists a unitary 

Fi *" 

where each Fi G Mm„. Then the map T defined 



U such that Y = UX. Write U ■ 

by T{A) = J^'.^^ FjAF* is a TPCP map sending XiX* to Bi for all i. 



: * 

Fr * 



n 
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Remark. To make the similarity to Theorem 13. II more apparent, note that both conditions 
in this theorem are equivalent to 

X*X = M o Y*Y for some correlation matrix M. (5) 

Indeed, if (|5]) holds, write M = C*C and C = [ci| . . . \ck]. Since Ma = 1, c, is a unit vector. 
Let Vi = Ci ® Hi and Y = [yi\ ■ ■ ■ \yk]- Then iji, i = 1, . . . ,k, are purifications of Bi and 
Y*Y = X*X, so we have the second condition in the theorem. The reverse implication is 
trivial. 

One definition for the fidelity between two states A and B is 

F{A,B) = WVaVbWi = suplltrv^v^V"! : V" is a contraction}. 

It is known that a necessary (but not in general sufficient) condition for the existence of a 
TPCP map sending Ai, . . . , A^ to Bi, . . . , B^ is that 

F{Bi, Bj) > F{Ai, Aj) for all l<i,j <k (6) 

(see [21 Lemma 1]). The corollary below allows us to deduce this fact immediately when 
the input states are pure (since F{xiX*,XjX*) = \x*Xj\). It also illustrates what missing 
information (namely, the partial isometrics Vi) is needed in conjunction with ([6]) to create 
a sufficient condition for the existence of a TPCP map. Unfortunately, it is still not very 
computationally efficient. 

Corollary 3.6. Suppose Xi, . . . , x^ G C" are unit vectors and Bi, . . . ,Bk G M^ are density 
matrices. Then there is a TPCP map T such that T{xiX*) = Bi for i = 1, . . . ,k if and only 
if there exist partial isometrics V-i G M^r such that 

\fBiViV*^i = B, and x*Xj = tr ^/B^^jVjV* for all i,j. (7) 

Proof. Suppose Vi, . . . ,Vk are partial isometrics satisfying ([7]). Let Yi = y/^Vi G M^^, write 
Yi = [yii\ ■ ■ ■ \yri], and define yi G C™ as in (g]). Then Bi = YiY* = Yl'j=iyjiyjv so yi is a 
purification of Bi. Since X*X = Y*Y for X = [xi| . . . \xk] and Y = [yi\ ■ ■ . \yk], the result 
follows by Theorem 13.51 



Conversely, by Theorem 13.51 we may assume there are purifications yi of Bi in the form 
of @ and X*X = Y*Y. Let Yi = [yul . . . \yri] G M^r- Since YiY* = Bi, there exist partial 
isometrics Vi G M^r such that Yi = y^Vi. But x*Xj = trF.*Y,-, so ([7]) holds. D 

4 General physical transformations on pure states 

Theorem 13.11 (quoted from [2j) gives a simple criterion for the existence of a TPCP map 
sending pure states xixl, . . . ,Xkxl to pure states yiyl, . . . ,ykyl. One might wonder how 
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to generalize this criterion to handle arbitrary interpolating completely positive (CP) maps. 
The remark in [2] after Theorem 7 seems to assert that there exists a CP map sending XiX* to 
yiy* if and only if X*X = MoY*Y for some positive semidefinite M (without any restriction 
on the diagonal entries of M). However, this condition is neither necessary nor sufficient. 

For example, let {61,62} be the standard basis for C^. Take xi = 61, X2 = {ci + 62) /V2 
and 1/1 = 61, 1/2 = 62. Then M o Y*Y = M o / is diagonal for any matrix M, but X*X has 
nonzero off-diagonal entries, so the condition is not satisfied. Nonetheless, there is a CP map 
sending XiX* to yiy*; let S G M2 be such that Sxi = yi. Then the CP map T{A) = SAS* 
works. 

On the other hand, let Xi = X2 = 61. Let yi = 61, 1/2 = 26i. Let M = (61 + 0.562)(6i + 
0.562)*. Then X*X = M o Y*Y is the matrix of all ones. But clearly there is no map T 
sending 616^ to both 616* and 46i6*. 

The following results consider interpolating CP maps and unital CP maps, generalizing 
Theorem 13. H and giving necessary and sufficient conditions in the same spirit as [2]. 

Theorem 4.1. Fix positive semi-definite rank-one matrices XiX* G M„ and y^y* G M^ for 
i = 1 . . .k. Let X = [xi|x2| . . . \xk] and Y = [yi\y2\ . . . \yk]- Then there exists a completely 
positive map T such that T{xiX*) = yiy* if and only if there exists a positive semi-definite 
matrix M e Mk with Ma = 1 such that kerX*X C ker M o {Y*Y). 

Proof. There exists a completely positive map T such that T{xiX*) = yiy* if and only if 

r 

3Fi,...,Fr e Mmn such that ^ FjXiX*F* = yiy* \/i = l,...,k 

i=i 

<^==^ 3Fi, . . . , Fj. G Mmn and unit vectors Ci, . . . , Cfc G C such that FjXi = Cijyi 

«^==^ 3Fj G Mm„,unit vectors Cj G C so that FjX = YTj where F^ is diagonal with (Tj)ii = Cij 

r 

<^=^ 3 diagonal F^ G M^ with /^P^F* = I^ such that rowspace YTj C rowspace X \/j 

i=i 

(equivalently, kerX C ker YTj Vj, or kerX*X C kei T*jY*YTj Vj) 
^^ kerX*X C kerM o Y*Y where {Mt)ij = (Tt)iiiTt)jj 

r 

and M = \^ Mf is a positive semi-definite matrix with Ma = 1 

n 

We will present a result on unital completely positive maps sending pure states to pure 
states as a corollary of the following more general result. Recall that for a rank r positive 
semi-definite matrix A G M„ with spectral decomposition A = XiUiul + ■ — h XrUrU*, where 
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{til, . . . ,Uk} C C" is an orthonorinal set of eigenvectors of A corresponding to the positive 
eigenvalues Ai, . . . , A^, the Moore-Penrose inverse A"*" of A has the spectral decomposition 
A+ = Ar^MiM* + ■ • • + X;^Uru;. 

Theorem 4.2. Fix rank-one matrices XiX* G M„ and yiy* G Mm for i = 1 . . .k. Fix a 
positive semi- definite matrix B G Mm- Let X = [a;i|a;2| . . . \xk] and Y = [yi\y2\ ■ ■ ■ Ivk]- Then 
there exists a completely positive linear map T : Mn — )■ Mm such that 

T{I) = B and T{xiX*) = y^y* for i = 1, . . . ,k, 

if and only if there exists a positive semi-definite matrix M G Mk with Ma = 1 such that 

(1) kerX*X C kerM o {Y*Y) and (2) Y[M o (X*X)+]r* < B, 

(with equality in (2) should X have rank n). Here X~*" denotes the Moore-Penrose inverse 
ofX. 

Proof. Note that the existence of a CP map T such that T{I) = B and T(xjX*) = yiy* is 
equivalent to the existence of Fi, . . . , F^ G Mmn satisfying 

r r 

(a) J2 FjX^x*F; = y,y* Wi and (6) J^ F,F* = B. 

Proof of Necessity: Condition (a) and the proof of Theorem 14.11 imply that FjX = YTj 
for some diagonal Tj G M^ with Yl^i=i^j^*j ~ ^k- Moreover condition (1) follows with the 
matrix M defined by Mij = YH=ii^t)ii0^t)jj- 

Let P denote the orthogonal projection XX'^, and let P-*" = /„ — P. Then FjP = 
FjXX+ = YTjX+, so 

r r r 

B = J2 F,F; = Y,{F,P + F,P^){PF* + P^F*) = J2 FjPPF* + F.P^F* 

j=i j=i j=i 

r r 

= Y^ YrjX+{x+yr*Y* + J] FjP^f* but x+(x+)* = (x*x)+ 

r 

= Y[M o {X*X)+]Y* + Y^ FjP^F* 

>Y[Mo{X*Xy]Y* 

with equality if P = /„, that is, if X has rank n. 

Proof of Sufficiency: Since M is positive semi-definite with Ma = 1, we can write M = 
C*C where C = [ci|c2| . . . \ck] G Mrk, and q is a unit vector for all i. If necessary, we may 
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append extra zero entries to the end of each q so that we may assume r > m. Define 
diagonal matrices Tt G M^ by {Tt)ii = Cu- Then 

r r r 

M o Y*Y = Y^ T*Y*YTj, M o (X*X)+ = J^ r^(X*X)+r*, and J] r^T* = 4. 

Condition (2) imphes 

B = Y[Mo {X*X)+]Y* + EE* for some E 

r r 

= J2 Yrj{X*X)+r*Y* + J2 GjP^G* 

where we may choose Gj G Mmk so that GjP-^G* is proportional to an eigenprojection for 
EE* with rank at most one. Note that P^ = if and only if X has rank n. 

Define Fj = YTjX+ + GjP^. Then 

r r 

Y^ EjF* = Y YTjX^{X^yT*Y* + GjP^G* + YTjX^P^G* + GjP^{X+yT*Y* 

r 

= Y YTjX^{X^yT*Y* + Gj-P^G* = B 
i=i 

since X"'"P-'- = X+(J — XX+) = 0, and the fourth term in the second sum is the adjoint of 
the third term. 

On the other hand 

FjX = YTj{X+X -/ + /) + GjP^X 

= -YTj{I - X+X) + YTj + Gj{I - XX+)X. 

But / — X^X is the orthogonal projection onto kerX; since condition (1) implies kerX C 
kerFr^- for all j, the first term must be 0. And (/ - XX+)X = X - XX+X = 0, so the 
third term vanishes. Thus FjX = YTj for all j, and 

r 

y^ FjXiX*F* = tjitj* for alH = 1, . . . , A; 
i=i 

as desired. D 

Corollary 4.3. Fix XiX* G M„ and yiy* G Mm for i = 1, . . . ,k. Write X = [a;i| . . . \xk\ and 

Y = [yi\ ■ ■ ■ \yk]- Then there exists a unital completely positive map T satisfying T{xiX*) = 
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Hiy* for all i = 1, . . . ,k if and only if there exists a positive semi-definite matrix M G Mk 
with Ma = 1 such that 

(1) kerX*X C ker M o {Y*Y) and (2) Y[M o {X*X)+]Y* < /„, 

(with equality in (2) should X have rank n). 

Proof. Take 5 = /„ in Theorem 14.21 D 

Corollary 4.4. Use the notation in Corollary \4.3[ There is a unital TPCP map sending 
Xix\., . . . , Xkxl. to yiyl, . . . , ykVl ^/ ^^'^ only if m = n and there exists a positive semi-definite 
matrix M G M^. with Ma = 1 such that 

(1) X*X = Mo {Y*Y) and (2) Y[M o (X*X)+]r* < /„, 

with equality in (2) should X have rank n. 

Note 

Reference [1| was brought to our attention by one of the referees. In it the authors 
independently obtain our condition (c) in Theorem 12. 2[ Moreover, they extend the result 
by allowing final states to have dimension greater than two, although it appears that our 
proof is self-contained, and uses more elementary methods. They also consider the problem 
of approximately mapping a set of initial states to a set of final states via CP maps. 
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